A new identity linking coefficients of a certain class of 
fiomogeneous polynomials and Gauss hypergeometric 

functions 



Philip W. Livermore, Glenn R. lerley 

Institute of Geophysics and Planetary Physics, Scripps Institution of Oceanography, UCSD, La Jolla, 

CA 92093 USA 



Abstract 

In investigating the properties of a certain class of homogeneous polynomials, we dis- 
covered an identity satisfied by their coefficients which involves simple 2F1 Gauss hyper- 
geometric functions. This result appears to be new and we supply a direct proof. The 
simplicity of the identity is suggestive of a deeper result. 
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The results that we report here arose when investigating an applied mathematics 
problem ostensibly unrelated to combinatorics: the self-generation of magnetic fields in a 
spherical geometry as applicable, for instance, to the Earth's fluid core. It happens that a 
certain set of homogeneous conditions apply to such magnetic fields, the so-called Taylor 
constraints, which arise in the physically interesting low-viscosity fast-rotation limit of 
the governing equations [H^ . These constraints take the form of the vanishing of certain 
integrals over surfaces of constant cylindrical radius of a quantity involving the magnetic 
field. We have recently proven that a necessary condition for any solution is the vanishing 
of a certain function of cylindrical radius on both inner and outer spherical boundaries of 
the core [2] . On adopting a truncated spatial discretisation based on spherical harmonics 
in solid angle and certain regular polynomials in radius, this amounts to the vanishing of 

Q{p.s)^^a'jk'Ijk{p,s) (1) 

regarded as a function of s at fixed values of p = 7/20 and p — 1. The coefficients ajk 
are real and 

Ijk{p,s)= I ' ' z^^iz-' + s^fdz. (2) 







In 11]), j and k belong to a prescribed finite set (whose size depends on the spectral 
truncation adopted in the numerical scheme) and s is the nondimensional cylindrical 
radius; the given values of p represent respectively the spherical radii of the inner and 
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outer boundaries of the Earth's fluid core. We note two points associated with any given 
(j, k). Firstly, the integrand in ((2]) is a homogeneous polynomial in (s, z) of degree 2{j+k); 
secondly, its limits are homogeneous of degree one. Thus Xjk{p,s) is homogeneous of 
degree 2{j + k) + 1 and can be written 

j-\-k 

It follows immediately that, for each choice of p, U]) becomes 

L 

Q{p,s)^^/p^Y.Mp)^'' (4) 
/=o 

where L is the maximum value oi j + k. 

We may therefore impose the required conditions, the vanishing of both Q(7/20,s) 
and Q(l,s), simply by setting each coefficient Ai appearing in ^ individually to zero. 
This procedure yields a set of constraints of size 2L. However, it happens that this set 
contains (in general) many degeneracies. Empirically, we found that the Ai, produced 
by choosing one of the two values of p as unity, are related by 

An{1) ^ An{p) + pM^N+m{p), N > max{j} (5) 

m— 1 

where p,m is proportional to a 2F1 Gauss hypergeometric, given below. In our Taylor- 
state problem, we will fix p — 7/20, but the above statement is true more generally. 
Thus, for any A'^ > max{j}, if the values of Ai appearing on the right hand side above 
are set to zero, then An{1) automatically vanishes and explicitly enforcing A]S[{1) = is 
unnecessary and introduces a degeneracy in the set of constraints. 

The source of the linear homogeneous condition ^ is immediately traced to the same 
condition on each Jjk individually, that is, 

1 „ .2 3+k~N 

Bn{1)=Bn{p) + ^^ Pn.{p) BN+,n{p), N>j (6) 

ni—1 

where Bi, defined in ([3]), is proportional to as follows from homogeneity. Note 

that the upper limit on the summation is the tightest possible since i3; = if / > fc + j. 
It is immediate that (O follows from © because the summand is independent of (j, k) . 

In this note we provide a proof of a generalization of (HI, and it is then a simple 
matter to use ([5|) to count the number of independent conditions related to enforcing 
Taylor's constraint. 

The particular structure of the integrand arises in the following manner. The mono- 
mial dependence of z stems from the alignment of the cylindrical integrals, defining 
Taylor's constraint, with the z-axis. The appearance of only through + comes 
about due to assumed C°° behaviour of the magnetic field: since both z and the square 
of the spherical radial distance, = + z^, are both C°°, any functional relation of 
the given form also inherits this property. Lastly, the form of the upper limit on in- 
tegration comes from the height of the intersection of a cylinder of cylindrical radius 

2 



s with a sphere of radius p. The elementary structure of the identity is suggestive of 
a deeper result, although the statement cannot be readily generalized, for instance, to 
odd exponents of z. It is possible that this result arises in other circumstances involving 
cylindrically symmetric constraints in fluid dynamics, for instance, the Taylor-Proudman 
condition in a rotating sphere . 

Immediately below we state the main result, but defer proof until two necessary 
Lemmas are stated and proved. 

Theorem 1. Let j and k be positive integers and let the coefficients Bi{p) be defined by 



■^0 (=0 
Then, for any positive integer N with k + j > N > j , the quantity 

Bn{p)-\ ^ turnip) BN+mip) (8) 

m— 1 

is independent of p where 

^_ 2r(m + l/2) /[l-m,l/2] ^ 

^"^^^ - V^TVTT)' ' I [l/2-m] 

Remarks 1. The statement in the theorem is considerably stronger than that of equa- 
tion ^ and relates expressions of the form ^ between any two values of p, not necessarily 
including unity. Additionally, since /ip = 2(1 — p^)~^, the theorem is equivalent to the 
p-independence of 

l-p2 ^+''^-^ 

9 ^ P'nBN+ m ■ 

?n— 

Although this is the most succinct form of the result, the statement in the theorem 
is easier to prove as iimip), for m > 1, are simply polynomials. The /i,„ have been 
arbitrarily normalised such that po = 1. 

Lemma 2. For any real values a, (3, the following are identities 

V 1 - fc - l/2\ /fc - l/2\ r(/3)r(n + /3+l/2) 

^„/3 + fcl -1/2 A -1/2 y r(/3+l/2)r(/3 + n+l)' ^ 

where (^) is a binomial coefficient (extended to non-integer values of its arguments). 

Proof. We rewrite the left hand side of the above two equations in a hypergeometric 
form, 

(-im«+i) f h. .+i] 



T{a - (3 - n + l)V{(3 + n + I) \[a~(3-n + l\ 

r(„ + i/2) f y ,,,, 



^r(n + l)^'' \ [-n+ 1/2,13+1]'' 



Using the theorems of Gauss and Saalschiitz [4|], we may evaluate the 2F1 and 3F2 
functions and the results follow easily. □ 

Proof of main result. We prove the theorem directly in two steps. First, we find a 
general expression for the coefficients Bi and second, we show the required relation holds 
between them. 

Step 1 By exploiting equation B6 in we can write 

The hypergeometric function appearing is simply a polynomial in / of degree 



2 



, 2 / 2\ _^ r(fc + i)r(fc + j-n + i/2) /g-"" 

M [l/2-j--fc]'" 'P ) ^^T{k-n+l)T{k + j + ll2) \p 



and Tjk{p, s) can be written 

IMP, s) = {p^ ~ oJ.^ltlL. - t "(^+^--"+^/^) .^» .^(-") 



2r(fc + j + 3/2) > r(fc-n + l) 



-a p 



where, aside from the leading prefactors, the binomial and summation appearing 
immediately above can be further written as 

V V r(fc + j-n + i/2)r(j + i)(-i)' 2(fc-„+,-.) ... 

1/2 

Up to the prefactor of (p^ — s^) , IjkXp, s) is simply a polynomial in of degree 
k + j. To express this in the form of ((3]), we may write I ^ n + i and re-order the 
sum over dummy indices / and n. However, care must be taken with the limits. 
Noting that i > and therefore n < I, coupled with the upper bound n < k, leads 
to n < min(fc,Z). Additionally, since i < j and therefore n > I — j , coupled with 
the lower bound n > 0, leads to n > max(0, 1 — j). Hence 



2{k+j-l) 



^ r(fc + i)r(j + i) ""^''^ (-i)'+"r(fc + j-n + i/2) 

^^P' 2r(fc + 7 + 3/2) ^ r(fc-n + i)r(7 + n-; + i)ra-n + i)^ 

(14) 

Shifting the dummy indices by writing m ~ n — (l ^ j) and defining T — k + j ~ I, 
this expression can be simplified to 

r(fc + l)r(j + l/2)p^(^+:'-0 "™ /T-l/2 + j-m\/ j 

m— inax(Oj— /) 

(15) 



Note that if m > T or m > j then respectively, the first or second binomial term 
vanishes. Hence, without loss of generality, we can replace the upper limit of the 
sum by J. If / > j then the lower limit is zero and a trivial reordering of the sum, 
setting m' = j — m, and using (O, leads to 

B,(p) = r(fc + J-1 + 1/2) r(fc + 1) r(, + 1/2) j-iy 

li I < j then the summation in (|15p has a nonzero lower limit and cannot be 
evaluated in such a simple closed form. 
Step 2 We now are in a position to prove that, if ^ > j then the Bi{p) satisfy ([H]). We shall 
show that 

' v^"^ R . , f r(fc + j-iv + i/2)r(fc + i)r(j + i/2) 

^^+™(/') - ^nk + j + l-N) T{k + J + 3/2) ^ ^ 

and it follows immediately that 

Bn{.P) + ^ ^ BN+m{p) ^J.rn{p) 

ni—l 

^ r(fc + j-A^ + i/2)r(fc + i)r(j + i/2) ( .(k+j-N) , _ 2^ v 
2V^r(fc + j + i-7V)r(fc + j + 3/2) P) P) 

^ r{k + j-N + 1/2) r(fc + 1) r(j + 1/2) 

2V^r(fc + j + l-7V)r(A: + j + 3/2) ' ^ ' 

independent of p. Note that ^ involves Bi with / > N. Coupled with the initial 
hypothesis N > j, this is consistent with I > j, required in the derivation of the 
closed form for Bi. 

It remains to show (|17p . By expanding the hypergeometric component of prnip), a 
polynomial of degree to — 1 in p, it follows easily that 

'^^^ 2 r(TO) Tjn + 1/2) T{m-n + 1/2) 
h'o'' r(m+l)r(m-n)r(n + l) ^ ' 

Hence 

''^'^o ,^ ,^ (-l)^T(fc+l)r(j + l/2) 



7r3/2r(fc + j+3/2) 

fe+^jv r(fc + J - TV - m + 1/2) r(m) r(n + 1/2) r(TO -n + 1/2) .(fc+^+n-iV-m) 

^ ^ r(/c + j + 1 - iv - TO)r(m + i)r(TO- n)r(n + 1) 

(19) 

This is a polynomial in of degree k + j — N — 1 . By introducing a new dummy 
variable b = k + j + n — N~m the double sum above can be rewritten as 



k+j~N k+j-N-1 

^ Bjv+rnip) fJ-mip) = 

m=l 6=0 



where 



_ (-i)j r(fc + 1) r(j + 1/2) T{k + j - N -b + 1/2) 

""^^ 7r3/2r(fc + j + 3/2)r(fc + j- Ar-5) ^ 

,n=k+^-N Y{k+j-N-m + 1/2) r(m) r(6 - fc - j + TV + m + 1/2) 

^ r(fc + 7 + 1 - AT - m) r(m + 1) r(5 - fc - j + iV + m + 1) ■ ^ ' 

The change in the lower limits on m arises since b — k— j + m + N = Ti>Om the 
original sum and hence m > k + j — N — b. Note that, since b<k+j — N — 1, 
m > 1, consistent with the original lower limit on m from (|19p . Defining new 
variables c = m — j3,j3 = k+ j — N — b, the summation can be written 

r(b + c+l/2)r(e+l/2) 1 ^ A /6 + c - l/2\ fc - l/2\ r(l/2)^ 
^ r(6 + c+l)r(c+l) c + /3"^lv -1/2 A-l/2y'c + /3 

and by using (llOp it is immediate that 

, r(fc + j-jv + i/2)r(fc + i)r(j + i/2) 

'''' ^ ^ ^T{k + ] + l-N)T{k + j + ?,/2) 
which is independent of 6, and p?)) follows. 
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